%calculates the local time from fday and longitude
%also calculates the LT at India from longitude
%latest date 31 OCT 2005 (B'day !)

function[lt] = champ_lt_new(fday,lon),

L = lon < 0;
lon(L) = lon(L)+360;

L = fday < 0;
ut(L) = (fday(L)-floor(fday(L)) + 1) * 24 ;
ut(~L) = (fday(~L)-floor(fday(~L)))* 24;
ut = ut';
 
lt = ut + lon/15;
L = lt > 24;
lt(L) = lt(L)-24;

    

